What drives the insulating state in ultrathin films of SrRuOs? 
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Using density functional calculations we have examined the evolution of the electronic structure 
of SrRuOs films grown on SrTiOs substrates as a function of film thickness. At the ultrathin limit 
of two monolayers (Ru02-terminated surface) the films are found to be at the brink of a spin-state 
transition which drives the system to an antiferromagnetic and insulating state. Increasing the film 
thickness to four monolayers, one finds the surprising result that two entirely different solutions 
coexist. An antiferromagnetic insulating solution coexists with a metallic solution corresponding 
to an antiferromagnetic surface and a ferromagnetic bulk. The electronic structure found at the 
ultrathin limit persists for thicker films and an unusual result is predicted. Thicker films are found 
to be metallic as expected for the bulk, but the magnetism does not directly evolve to the bulk 
ferromagnetic state. The surface remains antiferromagnetic while the bulk exhibits ferromagnetic 
ordering. 

PACS numbers: 75.70.-i,71.15.-m,71.30.+h,75.50.Ee 



There has been considerable efforts in the search for 
alternate technologies to the current Si-based technology 
[H . An avenue of research that has been intensively ex- 
plored is the use of transition-metal perovskite oxides. 
The interest in these materials stems from their multi- 
functional properties Several of these materials ex- 
hibit different types of magnetic ordering, ferroelectricity, 
charge and orbital ordering, all resulting from a strong 
interplay between spin, charge and orbital degrees of free- 
dom. A set of interesting members of this class of com- 
pounds are those in which the magnetic ordering can be 
manipulated by an electric field, or the ferroelectric or- 
dering can be manipulated by a magnetic field. Devices 
based on these materials have been envisioned and this 
has fueled the study of oxide electronics [sj . Metallic ox- 
ides form an integral part of the oxide-based device tech- 
nologies, as they would form the connecting electrode 
material. 

Among the perovskite oxides, SrRuOa is a promis- 
ing material device applications, which apart from be- 
ing a metal, although a poor one, also exhibits ferro- 
magnetism up to 160 K 4]. An additional advantage of 
SrRuOa is that correlation effects should be significantly 
weaker than in the 3d transition-metal oxides. However, 
when thin films of this material were epitaxially grown on 
SrTiOa, films less than four monolayers thick were found 
to be insulating 0, @] . While the earlier studies reported 
the loss of ferromagnetic ordering in ultrathin films 
recent studies have reported exchange bias effects which 
indicate the presence of antiferromagnetic regions Q . All 
in all, the exact mechanism driving the metal-insulator 
transition in SrRuOa ultrathin films is still unknown. 
Could the metal-insulator transition be due to electron 
correlations, or could it be due to structural effects or ex- 



traneous experimental conditions as recently suggested 
[7:]? The question of whether electron correlations are 
important in these systems has been investigated exper- 
imentally, although there is no consensus. Some exper- 
imental studies suggest that electron-correlation effects 
are important even in the 4(i compounds and cannot be 
ignored [1, 0, [l^. Other studies suggest that electron- 
correlation effects are much weaker and are not at all im- 
portant in this class of materials [ll!| . In order to address 
these issues, we performed first-principles density func- 
tional calculations of thin films of SrRuOa on a SrTiOs 
substrate, and examined the evolution of the electronic 
structure of the films with respect to the film thickness. 
We find that electron correlations at the level of GGA-|-?7 
must be included in order to explain the insulating state. 
The correlation effects serve the purpose of changing the 
subtle energy difference between the low-spin and the 
high-spin states. The transition into the high-spin state 
at the surface is accompanied by a strong Jahn- Teller- 
type distortion associated with a d'^ configuration with a 
partially filled eg states, which drives the system to an 
insulating state. The choice of U in GGA-\-U calculations 
is usually arbitrary. In order to reduce this arbitrariness 
we calculated U from first-principles using a technique 
recently developed by one of the authors [12[ . The mag- 
nitude of U determined drives the thin films insulating. 

SrRuOa is found to crystallize in an orthorhombic 
structure with four formula units per unit cell [l3i |. The 
Ru-O-Ru angles in the a-b plane are found to be around 
168 degrees. SrTiOa is found to occur in a cubic struc- 
ture. The GGA optimized lattice constant of ferromag- 
netic cubic SrRuOa is 3.99 A, while that of cubic SrTiOs 
is 3.94 A. Therefore, there exists a small mismatch be- 
tween cubic SrRuOs and SrTiOa. In order to simulate 
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the ultrathin films of SrRuOa grown on SrTiOa, we con- 
sider a 13 layer symmetric slab containing a central SrO 
layer, and alternating Ti02/SrO layers along the (001) 
direction as shown in Fig. 1. The SrRuOa/SrTiOs struc- 
tures were constructed to represent SrRuOa films de- 
posited on the SrTiOa substrate on a layer by layer mode. 
In every case we have included 15 A of vacuum, and we 
have checked that the results are converged with respect 
to the thickness of the vacuum layer. The in-plane lattice 
constant of the Ti02 layers as well as the Ru02 layers is 
set equal to the equilibrium lattice constant of SrTiOs. 
As the Ru-O-Ru angles in bulk SrRuOa deviate from 
180 degrees, we have allowed for a GdFeOs distortion of 
the RuOg octahedra and optimised the structure to find 
the minimum energy solution. Full optimization of the 
atomic positions is performed. In addition, a \/2 x \/2 
reconstruction of the (001) surface was considered, which 
allowed us to explore the possibility of antiferromagnetic 
ordering. 

The electronic structure of bulk as well as of thin films 
of SrRuOa were investigated using the plane-wave pseu- 
dopotential implementation of density functional theory 
and projector-augmented wave potentials in the VASP 
code [3[l^. The GGA approximation for the exchange- 
correlation functional was used. Electron-correlation ef- 
fects at the Ru sites were included through the GGA+U 
method within the Dudarev et al. formalism [l^ . where 
U represents the onsite Coulomb interaction strength and 
J represents the intra-atomic exchange. We used the val- 
ues oiU ~ 2.5 and J=0.4 eV for the Ru M states, which 
were determined for bulk SrRuOa by using a recently pro- 
posed scheme based on the random-phase approximation 
[12| . Although for the surface the U may be larger due to 
the reduced screening, we expect the qualitative aspects 
of the results to remain the same. A special k-point mesh 
of 6 X 6 X 6 was used for integrations over the the Bril- 
louin zone of the bulk, while a mesh of 4 x 4 x 6 was used 
for the slab calculations. A denser mesh of 8 x 8 x 4 was 
used in the density of states calculations for the slabs. A 
cutoff energy of 250 eV was used for the plane wave basis 
set. 

How does the electronic structure of bulk SrRuOs 
change in GGA+U?. We first examine the electronic 
structure of bulk SrRuOa as given by the GGA and 
GGA-I-J7 functionals. The comparison between Ru 4d 
partial density of states (PDOS) for U=0 (GGA) and 
U = 2.5 eV, J = 0.4 eV is shown in Fig. 2. The Fermi 
energy is used as reference by setting it to zero. We note 
that there is an increase in the exchange splitting for 
the finite U calculations. In order to examine the ori- 
gin of the enhanced exchange splitting we plot the up- 
and down-spin partial density of states in the insets of 
panels (a) and (b) in Fig. 2. The up-spin density of 
states which were contributing at the Fermi level now 
move deeper into the valence band and the system be- 
comes half-metallic as discussed earlier in the literature 



by Jeng et al. [l7|. In the down-spin channel we find 
that the effect of U is minimal. Hence the increased 
exchange splitting comes from the movement of up-spin 
states. Gross features of the partial density of states are 
otherwise similar and have been previously discussed in 
the literature [11]. For completeness, we repeat the dis- 
cussion here in order to contrast with the surface which 
shows a very unusual electronic structure. 

Ru in SrRuOa has a formal configuration. As the 
crystal-field splitting is quite large, the electrons occupy 
the t2g down-spin states after occupying the t2g up-spin 
states. The system appears to be at the brink of a tran- 
sition to a half-metallic state with a weak shoulder corre- 
sponding to the up-spin states at the Fermi level. It has 
been pointed out earlier that a larger value of U drives 
the system into a half- metallic state 17|. The moment 
on the Ru atom calculated by taking a sphere of radius 
0.9 A is 1 bohr magneton, significantly reduced from the 
ionic value. This difference can be explained by the size- 
able moments residing on the oxygen atoms. 

Having understood the electronic structure of bulk 
SrRuOa we proceed by examining the electronic struc- 
ture of one atomic layer of Ru02 grown on the SrTiOs 
substrate. Before discussing the results from the calcu- 
lations, we first speculate on what to expect. Strong 
crystal-field anisotropics at the surface are expected to 
result in a level ordering of doubly degenerate dyz and 
dxz levels followed by dxy levels. Thus if the system un- 
dergoes a transition to a nonmagnetic state the four elec- 
trons would go into the dyz, dxz levels. This would ex- 
plain the origin of the insulating ground state. With this 
simple model in mind, we examined various possible so- 
lutions. Along with the ferromagnetic configuration, we 
considered nonmagnetic as well as antiferromagnetic con- 
figurations. For U=0, the ferromagnetic configuration is 
the ground state, being more stable than the nonmag- 
netic configuration by 24 meV per Ru02 unit, whereas 
the antiferromagnetic solution is difficult to converge to. 
However, for U=2.5 and a J=0.4 cV on Ru, we find that 
the antiferromagnetic state is more stable by 49 meV per 
Ru02 unit. Therefore, in every case, magnetism survives 
at this ultrathin film limit. Most importantly, the transi- 
tion to the antiferromagnetic state is accompanied by a 
change in the electronic state from metal to insulator. 

What is the origin of the insulating state ? While it is 
clear from our results that the insulating state is linked 
to an antiferromagnetic order, the origin of the insulating 
state is puzzling, and we rely on an analysis of the pro- 
jected density of states to explain it. The orbital- and 
spin-projected density of states of Ru shown in Fig. 3, 
and the corresponding schematic energy-level diagram 
shown in Fig. 4. In the up-spin channel we find that 
the dxy, dyz and dxz states have a significant weight in 
an energy window near —5 eV as well as near eV, and 
the exact ordering of these t2c/-derived levels is not very 
clear. However, when we examine the contributions of 
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the down-spin counterparts, the ordering becomes clear. 
The states with dxz and dyz character are at lower ener- 
gies compared to the state with d^y character. 

The origin of this splitting of the t2g levels is clear when 
one examines the bond lengths in the RuOs unit at the 
surface. The in-plane bond lengths are equal to ~1.95- 
1.96 A, while the apical oxygen bond length is ^^2.13 A. 
Hence crystal field effects determine the ordering of the 
t2g states, with the dyz and dxz levels followed by d^y 
in each spin channel. This ordering is depicted schemati- 
cally in the left panel of Fig. 4. These levels interact with 
the corresponding symmetry levels on the oxygen atoms 
forming bonding and antibonding states. We illustrate 
this interaction explicitly for the up-spin t2g-derived lev- 
els in Fig. 3. As the charge transfer energy between the 
up-spin t2g and the oxygen atoms is small, the bonding 
states arising from this interaction have significant Ru d 
character. The hopping interactions however should be 
strongest for the in-plane dxy orbital. This results in an 
inversion of the ordering of the bonding levels generated 
as a result of the interaction. 

The strong crystal field that is present at the surface 
also induces a huge splitting of the e^-derived levels. The 
d^z^_j,2 states, which were unoccupied in the bulk, are 
now occupied, while the up-spin dx^-y^ state lies above 
the down-spin dxy, dyz and dxz states. Hence an unusual 
ordering of levels emerges from the broken symmetry at 
the surface. 

The next question we asked was what is the energetics 
involved in driving the unusual level ordering that one 
finds here. The generation of the large crystal distortion 
would cost enormous strain energy which we speculate 
must be compensated by the additional Hund's intra- 
atomic exchange present as a result of the high-spin state 
on Ru being favored. It should be noted that this crystal- 
field distorted structure exists above a particular value 
of U for the antiferromagnetic state alone. The non- 
magnetic and ferromagnetic solutions do not show such 
distortion at the two monolayers limit. 

Experimentally it is found that the insulating state 
exists for four atomic layers, so it is important to ask 
whether the insulating state survives when we have an 
additional Ru02 layer. This has direct correspondence 
with experiment as it corresponds to four atomic layers 
on the exposed Ti02 surface. The Ti02 surface is usually 
the one exposed in the experiments on which the SrRuOs 
films are grown [l^. The Ru d-projected PDOS for the 
atoms at the surface layer as well as at the sub-surface 
layer are shown in Fig. 5. The system is still found to 
be insulating. Interestingly the surface Ru atoms show a 
similar structural distortion as found in the single Ru02 
layer case. However, the inner Ru02 layer does not show 
the same distortion but still is insulating. Thus it is the 
antiferromagnetic correlations which set in as a result of 
the crystal distortion of the surface layer which drives 
the sub-surface Ru02 layer insulating. 



In order to probe the thickness at which we have an 
insulator-metal transition, we investigated the effects of 
additional Ru02 layers on the electronic structure of the 
films. Surprisingly, we find that the antiferromagnetic in- 
sulating solution is always lower than the ferromagnetic 
metallic solution for even 8 monolayers. This is far be- 
yond the thickness for which a metal-insulator transition 
has been observed in experiments. One might suggest 
that U—2.5 eV used in our calculations is too large for 
the ruthenates. However, we note that the exchange bias 
effects that has been recently observed in these systems 
offer us some hints to the solution of this puzzle. Could 
the experimental results be explained by an antiferromag- 
netic surface and a ferromagnetic bulk (AFMS-FMB) ? 
By examining this further we found that the antiferro- 
magnetic solution is actually degenerate with the AFMS- 
FMB state at the four monolayer limit. The antiferro- 
magnetic solution as we pointed out earlier is insulating, 
while the AFMS-FMB solution is metallic. Thus, our re- 
sults indicate that two electronically different states co- 
exist, and that disorder would possibly pin one solution 
in one region. For thicker films of six or more monolayers 
we find the AFMS solution to be the most stable. The 
AFMS-FMB solution we find is the favored solution far 
away from the ultrathin limit, even for SrRuOa surfaces. 

Overall, thin films of SrRuOa of 4 monolayer thickness 
or less are found to be insulating. Examining the elec- 
tronic structure of ultrathin films with just one Ru02 
layer, we find that the Ru atom at the surface undergoes 
a low-spin to high-spin transition. This rare occurrence 
of a high-spin state in a 4(i oxide results in an antiferro- 
magnetic state being stabilized as the ground state. The 
antiferromagnetic ordering drives the system insulating. 
In the limit of 4 monolayers, the antiferromagnetic insu- 
lating solution coexists with the solution corresponding 
to an antiferromagnetic surface and a ferromagnetic bulk 
(AFMS-FMB). For six-monolayers thick films, we find 
antiferromagnetic surface and the bulk ferromagnetic so- 
lution, which is metallic, to have lower energy. In other 
words, we find that SrRuOa surfaces are antiferromag- 
netic, while the bulk is ferromagnetic. 

Hence the electronic structure of systems in the ultra- 
thin limit are very different from the bulk. Crystal-field 
anisotropics largely determine the properties of ultrathin 
films. In the case of SrRuOa we show that a metal to 
insulator transition takes place as a function of thick- 
ness in addition to a ferromagnet to an antiferromagnet 
transition. 

PM thanks the Department of Science and Technology 
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FIG. 1: A schematic diagram of the slab geometry that 

was considered. The substrate was fixed at 11 monolayers 
on which the SrRuOs films were grown. The present figure 
corresponds to 2 monolayers of SrRuOs. 
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FIG. 2: The up (solid line) and down (dashed line) spin Ru d 
partial density of states for bulk SrRuOs in the ferromagnetic 
state. The zero of energy corresponds to the Fermi energy. 
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FIG. 3: The up (solid line) and down (dashed line) spin or- 
bital projected Ru d partial density of states for the surface 
Ru atom for 2 monolayers of SrRuOa grown on SrTiOa cal- 
culated for the antiferromagnetic state. The zero of energy 
corresponds to the top of the valence band. 



Energy level diagram of RuOg terminated surface 
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FIG. 4: The up (solid line) and down-spin (dashed line) 

crystal field split Ru d levels for the surface atom are shown 
schematically in the left panel. The interaction of these levels 
with the O p orbitals (in the right panel) and the ensuing 
ordering of the up-spin t2g-derived levels are shown in the 
central panel. 
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FIG. 5: The up (solid line) and down (dashed line) spin Ru 
d partial density of states for 4 monolayers of SrRuOa grown 
on SrTiOa for (a) sub-surface Ru atom as well as for (b) the 
surface Ru atom are shown for the antiferromagnetic state. 
The zero of energy corresponds to the top of the valence band. 



